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Abstract 

The volume flow rates of Eyring-fluids inside the wavy-rough nanoannulus were obtained 
analytically (up to the second order) by using the verified model and boundary perturbation 
method. Our results show that the wavy-roughness could enhance the flow rate especially for 
smaller forcing due to the larger surface-to-volume ratio and slip-velocity effect. Meanwhile, 
the phase shift between the outer and inner walls of nanoannuli could tune the transport 
of Eyring-fluids either forward or backward when the wavy-roughness of a nanoannulus is 
larger enough. Our results could be applied to the flow control in nanofluidics as well as 
biofluidics. 

PACS : 83.60.Rs, 83.50.Lh 

1 Introduction 

Solutions of linear chain structures exhibit interesting rheological properties such as the shear 
thinning and suppression of turbulent flow, which is related to flow induced changes in chain con- 
formation and orientation. These flow properties are of great relevance in technical applications 
as thickeners, drag reducers, and flow improvers, as well as in the production of fiber-reinforced 
materials. A unique property of long chainlike molecules is the formation of entanglement net- 
works. Already at low concentrations chain molecules start to overlap and entangle to form a 
transient network. Chainlike molecules subjected to a viscous shear gradient will orient in the 
flow, the instantaneous angular velocity being a function of the orientation relative to the local 
streamlines [I]. 

Meanwhile, researchers have been interested in the question of how some material responds to 
an external mechanical load [2]. External loads cause liquids to flow, in Newtonian or various 
types of non-Newtonian flows. Glassy materials, composed of polymers, metals, or ceramics, can 
deform under mechanical loads, and the nature of the response to loads often dictates the choice 
of material in various industrial applications. In biological systems, the response of proteins to 
external loads governs aspects of cell adhesion and muscle function [3]. The nature of all of 
these responses depends on both the temperature and loading rate. As described by Eyring [3], 
mechanical loading lowers energy barriers, thus facilitating progress over the barrier by random 
thermal fluctuations. The Eyring model approximates the loading dependence of the barrier 
height as linear. The Eyring model, with this linear barrier height dependence on load, has been 
used over a large fraction of the last century to describe the response of a wide range of systems 



[5] and underlies modern approaches to biophysical rupture processes [6], sheared glasses [7], etc. 
To the best knowledge of the author, the simplest model that makes a prediction for the rate and 
temperature dependence of shear yielding is the rate-state Eyring model of stress-biased thermal 
activation [3,7]. Structural rearrangement is associated with a single energy barrier E that is 
lowered or raised linearly by an applied stress a. In glasses, the transition rates are negligible 
at zero stress. Thus, at finite stress one needs to consider only the rate ii^ of transitions in the 
direction aided by stress. 

The linear dependence will always correctly describe small changes in the barrier height, since 
it is simply the first term in the Taylor expansion of the barrier height as a function of load. 
It is thus appropriate when the barrier height changes only slightly before the system escapes 
the local energy minimum. This situation occurs at higher temperatures; for example, New- 
tonian flow is obtained in the Eyring model in the limit where the system experiences only 
small changes in the barrier height before thermally escaping the energy minimum. As the 
temperature decreases, larger changes in the barrier height occur before the system escapes the 
energy minimum (giving rise to, for example, non-Newtonian flow). In this regime, the linear 
dependence is not necessarily appropriate, and can lead to inaccurate modeling. For example, 
Li and Makarov [8] have shown that there is a nonlinear barrier height dependence in stretched 
proteins, and that the assumption of a linear dependence in the analysis of experimental results 
leads to inaccurate conclusions. To be precise, at low shear rates (7 < 7c), the system behaves 
as a power law shear-thinning material while, at high shear rates, the stress varies afflnely with 
the shear rate. These two regimes correspond to two stable branches of stationary states, for 
which data obtained by imposing either a or 7 exactly superpose. The transition from the 
lower branch to the higher branch occurs through a stable hysteretic loop in a stress-controlled 
experiment [9]. 

Note also that one prominent difference between the fluid motions in nanodomain and those in 
macrodomain is the strong fluid-wall interactions observed in nanoconduits. For example, as the 
nanoconduit size decreases, the surface-to-volume ratio increases. Therefore, various properties 
of the walls, such as surface roughness, greatly affect the fluid motions in nanoconduits. 
In this short paper, we adopt the verifled Eyring model [3-4] to study the transport of shear- 
thinning fluids within corrugated nanoannuli. To obtain the law of shear-thinning fluids for 
explaining the too rapid annealing at the earliest time, because the relaxation at the beginning 
was steeper than could be explained by the bimolecular law, a hyperbolic sine law between the 
shear (strain) rate : and (large) shear stress : r was proposed and the close agreement with 
experimental data was obtained. This model has sound physical foundation from the thermal 
activation process [3-4] (Eyring [3] already considered a kind of (quantum) tunneling which re- 
lates to the matter rearranging by surmounting a potential energy barrier). With this model 
we can associate the (shear-thinning) fluid with the momentum transfer between neighboring 
atomic clusters on the microscopic scale and reveals the atomic interaction in the relaxation of 
flow with dissipation (the momentum transfer depends on the activation shear volume, which 
is associated with the center distance between atoms and is proportional to fc^T/ro (T is tem- 
perature in Kelvin, and tq a constant with the dimension of stress). Thus, this model could be 



applied to study transport of shear-thinning fluids in nanodomain [10]. 

To consider the more realistic but complicated boundary conditions in the walls of nanoan- 
nulus, however, we will use the boundary perturbation technique [11] to handle the presumed 
wavy-roughness along the walls of nanoannuli. The relevant boundary conditions along the 
wavy-rough surfaces will be prescribed below. 

2 Physical Formulations 

We shall consider a steady transport of the (shear-thinning) fluids in a wavy-rough nanoannulus 
of r2 (mean-averaged outer radius) with the outer wall being a fixed wavy-rough surface : 
r = r2 + esin(k9 + (3) and ri (mean-averaged inner radius) with the inner wall being a fixed 
wavy-rough surface : r = ri + €sm{k9), where e is the amplitude of the (wavy) roughness, /3 
is the phase shift between two walls, and the roughness wave number : k = iir/L. Firstly, this 
fluid [3-4,10] can be expressed as 7 = 70 sinh(r/To), where 7 is the shear rate, r is the shear 
stress, and 70 is a function of temperature with the dimension of the shear rate. In fact, the 
force balance gives the shear stress at a radius r as r = —{rdp/dz)/2. dp/dz is the pressure 
gradient along the flow (or tube-axis : z-axis) direction. 

Introducing x = —{f2/'2To)dp/dz then we have 7 = 70 sinh(xr/r2). As 7 = —du/dr {u is the 
velocity of the fluid flow in the longitudinal (2;-)direction of the nanoannulus), after integration, 
we obtain 

u = Us-\ coshy — cosh( — ), (1) 

X r2 

here, Ug is the velocity over the (inner or outer) surface of the nanoannulus, which is determined 
by the boundary condition. We noticed that Thompson and Troian [12] proposed a general 
boundary condition for transport over a solid surface as 

= l07(1 - ;^)-^/^ (2) 
7c 

where Au is the velocity jump over the solid surface, L° is a constant slip length, % is the critical 
shear rate at which the slip length diverges. The value of 7c is a function of the corrugation of 
interfacial energy. 

With the boundary condition from Thompson and Troian [12], we can derive the velocity fields 

and volume flow rates along the wavy-rough nanoannulus below using the verified boundary 
perturbation technique [11]. The wavy boundaries are prescribed as r = ri + esm{k6) and 
r = r2 + e sm{k9 + /?) and the presumed steady transport is along the z-direction (nanoannulus- 
axis direction). 

Along the outer boundary (the same treatment below could also be applied to the inner bound- 
ary), we have 7 = (du)/{dn)\Q^ surface- Here, n means the normal. Let u be expanded in e : 
u = uq + eui + €^U2 + • • •, and on the boundary, we expand u{ro + edr, 6{= Oq)) into 

dr^ 

u(r,e)\(^ro+€dr,eo) = u{ro,0) + e[dr Ur{ro,9)] + e^[—Urriro,9)] H = 

{usiip + — [coshx - cosh(— )]}|on surface- ^0 = ri,r2; (3) 

X ^2 



where 



""■siip I on surface ~ -^5 • ) ^ ]} Ion surface' (^) 

Now, on the outer wall (cf. [11]) 

^ _ du _ V(r - r2 - esin{k9 + /?)) 
dn |V(r — r2 — esm(A:6' + /?))! 

Considering ~ ri,r2 ^ e case, we also presume sinhx ^ ic/io- With equations (1) 
and (5), using the definition of 7, we can derive the velocity field (u) up to the second or- 
der : u{r,9)= — (r27o/x){cosh(xr/r2)— coshx [1 + e^X^ sin^(/c^ + /?)/(2r2)]-|-exsinhx sin(A;0 + 
siip\r=r2+esm{kd+f3)- The key point is to firstly obtain the slip velocity along the bound- 
aries or surfaces. After lengthy mathematical manipulations, we obtain the velocity fields (up 
to the second order) and then we can integrate them with respect to the cross-section to get the 
volume fiow rate {Q, also up to the second order here) : Q = /q^ Iri+esin{k9)^^^ ^)'^drd6 = 
Qsiip + ^Qpo + Qp2- In feet, the approximate (up to the second order) net volume fiow rate 
reads Q = Qout — Qin which is the flow within the outer (larger) wall : Qout without the contri- 
butions from the flow within the inner (snaller) wall Qm- 



3 Results and Discussions 

We shall demonstrate our results below. The wave number of roughness is fixed to be 10 
(presumed to be the same for both walls of the nanoannulus) for all figures here. Firstly, 
there is an enhanced fiow rate once the wavy-roughness is increasing [10]. This enhancement is 
rather significant especially when the forcing (along the annulus-axis direction : xl'^'i) is absent 
or zero (purely slip flow). As the annulus size decreases, the surface-to- volume ratio increases. 
Therefore, surface roughness along both walls together with the slip- velocity boundary condition, 
greatly affect the fluid motions between the corrugated nanotubes. 

Not that for an easy comparison, we select the parameters to be r2 = 2,ri/r2 = 0.5,Lg/ri = 
1; 7o/7c = 0.1. For rather weakly corrugations : e = 0.06ri the fiow rate is monotonically 
decreasing for around 5. However, once we increase to be around 7, the fiow rate will 
be monotonically increasing. Q (the net volume flow rate) will firstly decrease to a minimum 
and then keep increasing monotonically. This behavior is almost the same for /? (phase shift) 
being equal to 7r/4, 7r/2 and vr (however, the effect of phase shift is minor for very-small wavy- 
roughness [10]). It seems to us that, however, for very-small applied forcing there is a barrier 
manifested by the shearing (of fluids along wavy-rough surfaces) which leads to a minimum fiow 
rate. As the forcing is large enough and the barrier can be overcome then the transport (of 
fluids) keeps increasing. The latter is something like a tunneling process [3-4]! 
To examine the realistic effects of phase shift (/?) between the outer-wall wavy-roughness and the 
inner- wall one, we flx the outer- wall and the inner- wall radii (to be 1 and 0.6 nm, respectively). 
The small amplitude of wavy-roughness of both walls is also fixed to be O.OGri. We consider two 
cases : /? = 7r/4, vr for the same 70 = 10000.0 (s^^) (cf. [H]) with 7o/7c = 0.1. The results are 



illustrated in Fig. 2. We observe that for /3 = 7r/4 the net volume flow rate Q = Qout — Qin is 
decreasing as the forcing increases (starting from zero but within a small range). The trend for 
/? = TT, however, reverses! Q for /3 = tt is monotonically increasing and positive for larger forcing. 
Here, the interesting observation is for small applied forcing (4.5 < x/r2 < 7.5, r2 is fixed to be 
unity), once j3 = n/A, the flow moves backward (along the annulus-axis direction, even though 
the flow still moves forward once x/f2 > 8 [14] as also evidenced in Fig. 1 since there is a 
minimum flow rate for larger wavy-roughness for certain under selected geometry). This 
result thus could be applied to the flow control in nanofluidics. 

In brief summary, we have theoretically obtained the volume flow rates (up to the second order) 
of Eyring-fluids inside the wavy-rough annular nanotubes by using the verified fiuid model [3- 
4,10] and boundary perturbation method [11]. Our results show that the wavy-roughness could 
tune the flow rate especially for smaller forcing due to the larger surface-to-volume ratio and 
slip-velocity effect. Meanwhile, the phase shift between the outer and inner walls of nanoannuli 
could tune the transport of shear-thinning fluids either forward or backward when the wavy- 
roughness is larger enough as illustrated in Figure 2 here. Our results could be applied to the 
flow control in nanofluidics [15] as well as biofluidics [1]. 
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Fig. 1 Calculated net volume flow rates (Q) w.r.t. x/f2 (forcing (along the z-axis 
direction) per unit volume and referenced shear stress). The mean outer radius is 
r2 = 2 (nm), and the mean inner one is ri = 1 (nm). The slip length reads = ti. 
e (=0.06 ri here) is the amplitude of wavy roughness. The wave number of 
roughness (k) is 10 here. We demonstrate the effects of wavy-roughness via its phase shift 
/3 (=7r/4, 7r/2, it here) between the outer and inner walls of a nanoannulus. For smaller /?, 
Q is smaller w.r.t. the same x/f2 (forcing). Meanwhile there is a minimum flow rate 
(corresponds to a barrier). The unit of volume flow rate is mP/s. 
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Fig. 2 Calculated volume flow rates (Q) w.r.t. xlr2 (forcing (along the z-axis 

direction) per unit volume and referenced shear stress) for different phase shifts (/3s). 

The mean outer radius r2 = 1 nm, and the mean inner radius r\ = 0.6 nm. 

The amplitude of wavy-roughness e (=0.06 r\ here) is enlarged here and the wave 

number of roughness (A:) is 10 here. We consider two phase shifts (/? = 7r/4, vr) and 

the effect is significant for the range of here. The unit of volume flow rate is / s. 

Q for /3 = 7r/4 is monotonically decreasing (up to dimensionless x/'^2 > 8 as r2 is unity) 

and there is backward transport for this case. The trend for (3 = tt reverses. 

This result could be applied to the flow control in nanofluidics. 



